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ABSTRACT 



Context. Multiwavelength observations of supernova remnants can be explained within the framework of diffusive shock acceleration 

theory, which allows effective conversion of the explosion energy into cosmic rays. Although the models of nonlinear shocks describe 

reasonably well the nonthermal component of emission, certain issues, including the heating of the thermal electron plasma and the 

related X-ray emission, still remain open. 

Aims. Study of electron heating in supernova remnants in the Sedov phase due to Coulomb exchange between protons and electrons. 

Methods. Numerical solution of the equations of the Chevalier model for supernova remnant evolution, coupled with Coulomb 

scattering heating of the electrons. 

Results. The electron temperature and the X-ray thermal Bremsstrahlung emission from supernova remnants have been calculated 

as functions of the relevant parameters. Since only the Coulomb mechanism was considered for electron heating, the values obtained 

for the electron temperatures should be treated as lower limits. Results from this work can be useful to constrain model parameters 

for observed SNRs. 



1. Introduction 

Multiwavelength observations of supernova remnants (SNRs) 
allow one, in certain cases, to determine the temperatures Tg 
and Tn of both electrons and protons behind the shock, (see 
iRak owskr 2005, for a recent review). The ratio T^/Tp of these 
temperatures varies significantly, from < 0.07 for SN 1006 
Caming et al. 1996, Ghavamian et al. 2002) to 0.67 - 1 in the 
Cygnus Loop (.Ghavamian et al..,2001,) . From a theoretical point 
of view there are two limits for this ratio. The lower limit, given 
by the electron to proton mass ratio T^/Tp — nie/mp, is the so 
called total non-equilibrium case, which corresponds to equal 
thermal velocities for electrons and protons; the upper limit, 
Te = Tp is the total equilibrium case, which corresponds to 
eq ual energies of el e ctron and p roton plasmas. As was shown 
bv lGhavamianetal.1 (1200 ll 12002 ) for the SNRs Cygnus Loop, 
RCW 86, Tyco and SN 1006, observed values lie between these 
two limits. A compilation of current electron-ion eq uilibration 
measurements at SNRs shocks has been presented in lRakowskil 
(l2005h . where an inverse relationship between the shock veloc- 
ity and the level of temperatures equilibration has been claimed. 
Namely, older SNRs, characterized by a lower value of the shock 
velocity, exhibit a T^/Tp ratio closer to, and in some cases con- 
sistent with 1 . 

It is important to stress that, even in the case of total non- 
equilibrium at the shock {T^/Tp = m^/mp), electrons are 
unavoidably heated sufficiently far downstream of the shock 
throu gh Coulomb collisions with hot thermal protons dSpitzeiJ 
11998 ). This constitutes a sort of minimal heating scenario for 
electrons, which will be studied in details in this paper. The tem- 
per ature of this electrons is measured directly fr om X-ray spectra 
(e.g lDecourcheUe et allllooitlVink et alJl2003h . 



However, the electron heating through Coulomb interactions 
has been claimed to be too slow to explain the level of equilibra- 
tion measured from optic al lines observ ations in the vicinity of 
the shock for some SNRs (lLamingll2000l) and other mechanisms 
based on the generation of plasma w aves are belie ved to pro- 
ceed faster than Coulomb heating (e.g. Laming'200r). Several of 
such collisionless mechanis ms have b een inves tigated, a mongs t 
others, by ICa rgill & Papa dopoulosI ( 1988): Laming! (l200lh : 
iGhavamianeTalj (i2007i) : iRakowski et alj (i2008i) . 

The efficient acceleration of cosmic rays at SNR shocks also 
plays a role, by reducing the total energy available for shock 
heating and cons equently supp ressin g the proton temperature 
(Blasi et al."2005'; 'Ellison et al."2007'; 'Vink^OOSt iHelder et al.l 
2009; Drury et al. 2009; Patnaude et al. 2009). 

The aim of this work is to present a simple model for elec- 
tron Coulomb heating at SNR shocks, that could help to estimate 
the need for, and magnitude of, heating processes adjunctive to 
the Coulomb mechanism. Our work is based on Chevalier's self- 
similar model which describes the evolution of SNRs i n presence 
of eff ective cosmic ray acceleration at the SNR shock dChevalieiJ 
Il983h . We coupled Chevalier's model with the equations for 
electron heating by Coulomb exchange. The model lacks some 
important features; it does not include magnetic fields explicitly 
and the shock compression ratio is limited to a maximum value 
of 7 (since it does not include possibility for cosmic rays to es- 
cape from the shock), whereas in numerical simulations this ratio 
in some cases could reach much higher values (see for example 
Ellison et al. 2005). However it does include a minimal set of the 
key physical processes and provides, we feel, a useful reference 
calculation for more complicated estimates and models. 



D. Malyshev et al.: A simple model for electron plasma heating in supernova remnants. (RN) 



The paper is organized as follows. In section |2] we briefly 
summarize the Chevalier model, introduce the equation for heat- 
ing, and discuss projection effects in relating the 3D model to 
2D observations. In section [3] we summarize the results of the 
numerical solution of the equations of our model, and in section 
|4]we apply the results of our model to the case of SN 1006. We 
conclude in Sec. ID 

2. SNR dynamics 

2.1. Chevalier model 

A formal approach to the hydrodynamics of a SNR in adiabatic 
expansion ph ase includ i ng c osmic rays (CRs) acceleration was 
developed bv lChevaUed (Il983.) . 

In his work Chevalier generalized the self-similar model of 
Sedov, which describes a point explosion in a uniform medium 
with negligible pressure. Chevalier's model describes the inte- 
rior of a spherically symmetric SNR by the two-fluid (relativis- 
tic cosmic rays and non-relativistic thermal gas) hydrodynamical 
equations without diffusion (including diffusion would in gen- 
eral break the self-similarity). In terms of the radial coordinate 
< r < rshock{t) the equations are: 



dp dp dv 2pv 
ot or or r 



(1) 



dv dv 
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Here p and v stand for density and velocity of the gas, Pc and 
Pg pressures of the CRs and the thermal gas correspondingly. 
This system of equations describes mass and momentum con- 
servation (first two equations) and entropy conservation for gas 
and cosmic rays (last two equations). 

The boundary conditions come from applying generalized 
(i.e. cosmic rays included) jump conditions at the shock, which 
come from mass, momentum a nd energy conservation for par- 
ticles that cross the shock, see IChevalied (Il983h for details. To 
preserve self-similarity, the shock compression ratio, s, and the 
proportion of the total pressure going to cosmic rays w (CR ac- 
celeration efficiency) must be taken as constants giving: 



p{r shock) = spo 

v{r shock) =1 ) r shock 



(2) 



Pc{r shock) =w [1 ] po 



shock 



Poirshock) = (1 - w) ( 1 - - ) Por'ihock 



In fact s and w should be taken from a nonlinear shock model, 
but at the level of this model it is enough to treat them as ad- 
justable parameters. 

In standard fashion the system of partial differential equa- 
tions ([T]i can now be reduced to a system of ordinary differential 
equations by transforming to the similarity variable ^ ex rt^^^'^, 
with the constant of proportionality chosen such such that ^ ~ 
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utions of Chevalier model (see iChevi 
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983|) 



density n, velocity v and pressure Pq inside the SNR. Different 
curves coiTesponds to different values of CR acceleration effi- 
ciency w. Here ^ = r/r shock 



f/rshock{t) (see IChevalierl (Il983h for the constant's numerical 
value). 

Introducing scaled structure functions G, U and Z we obtain 



p(r,t)-poG(0; 
v{r,t) = \-U{(,): 
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for the density, velocity, and pressure, respectively, and for the 
compression ratio at the shock we get s = q(IZi_1 = '^"^j^ , 
where 7^ =3^ 

In the following we will work in the coordinate system de- 
fined by the variables ^ and t, so an element of fluid intersects the 
shock (^ = 1) and then moves with velocity ^ = |^(t^(C) ^ 1) 
towards the center (^ == 0) of the SNR. 

The solution of these equations for different w are shown in 
FiglT] where n = p/nip is the particle number density. 



2.2. Electron heating 

We now consider Coulomb electron interactions with hot pro- 
tons in the SNR model described in the previous subsection. 
In the coordinate system that moves with the plasma, the equa- 
tion th at describes electron heating, can be found e.g. in Spitzed 
(Il998l) . In the case of adiabatic expansion this equation reads: 



dt 



T —T 



7-l y dn. 
n dt ' 



(4) 



te, ^Ti'/'lnA/Tio; 



2.5 



1.5 



0.5 



Predicted flux 
Thermal bremsstralung spectrum 




0.2 0.4 0.6 O.i 



1 
E.keV 



1.2 1.4 1.6 



Fig. 2. The points show the predicted spectra from area dS at 
projected distance x — 0.7 from the SNR center for a SNR 
age t = 500 yrs, explosion energy E = 10^^ erg, gas density 
n = 0.1 cm^"^, CR acceleration efficiency w = 0. The solid line 
represents a single temperature TB spectra. 



where hi A is the Coulomb logarithm. 

The first term in Eq.|4]corresponds to the heating of electrons 
due to Coulomb interactions with the protons of temperature Tp 
and the second one corresponds to the cooling of the electron 
gas due to the adiabatic expansion of the SNR. This equation 
has to be coupled with the equation that connects Tp and Tg with 
macroscopic parameters such as the gas pressure and density and 
with the equation that allows us to change coordinate system 
from the system moving with the gas to the system connected 
with the center of the explosion: 

n{0 {Te{^,t) +Tp{^,t)) ^ PaiU); (5) 

P{£.) = mpn{^) 



Together with the boundary condition Te{£, = l,t)/Tp{£_ = 
1, i) = Trie /nip this forms a complete set of equations for de- 
termining the electron temperature Te{^, t) inside a SNR. 

2.3. Projection effects 

By solving Eq. |4] we derive the electron temperature Te{£,,t) 
in the radial coordinate ^. In order to compare this profile 
with observations we have to take into account projection ef- 
fects. To do this we use the following formula for the thermal 
bremsstr ahlung (TB) emission of the volume element dV at dis- 
tance D dRvbicki & Lightmanlll986l) : 
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where g// - is the thermal Gaunt factor, g// w 1, and the elec- 
tron temperature Tg here is a function of the distance ^ to the 
SNR center. 
The integration of this expression along the line of sight gives the 
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Fig. 3. Electron temperature determined from the peak of the x- 
ray emission (see Fig|2]l as a function of the (projected) distance 
from the center of the SNR for different explosion energies E, 
ages t and CR acceleration efficiencies w. The gas density is 
n = 0.1 cm^'^ for all curves. Triangles indicate the temperature 
obtained by integrating the X-ray emission from the whole SNR. 



flux from an area dS at a projected distance x = '"/'^s/i from the 
SNR center. The red points in Fig. |2]represent the expected flux 
as a function of energy for x = 0.7. As it can be easily seen this 
spectrum is not exactly thermal, being the superposition of dif- 
ferent thermal spectra characterized by different temperatures. 
However, a thermal bremsstrahlung spectrum gives a very good 
approximation to the predicted points. This is shown in Fig. |2] 
where the solid line represents the thermal bremsstrahlung emis- 
sion for an effective temperature defined as the peak of the ex- 
pected x-ray emission. The dependence of this effective temper- 
ature from X is shown in Fig. [3] 

This figure shows that, for a broad range of parameters such 
as explosion energy, SNR age and CR acceleration efficiency, 
such an effective electron temperature is almost constant in the 
inner (x < 0.9) region of the SNR. So, it is reasonable to study 
the dependence of Tg on other parameters not for all the values 
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Fig. 4. Predicted total thermal flux profile for a SNR of age 
1000 yr, explosion energy 10^^ erg expanding in an ambient 
medium with density 0.1 cm"'^. The effects of CR acceleration 
has been neglected (w — 0). Solid red, dashed green and dot- 
dashed blue lines refer to the flux at energies equal to 0.5, 1 and 
2 keV, respectively. 



of X, but only for some fixed x < 0.9. In the following we will 
adopt X — 0.7. Finally, the triangles in Fig.|3]show the tempera- 
ture one would observe by integrating the whole X-ray emission 
from the SNR. The discrepancy between the triangles and the 
curves is always below 5%. 

In Fig. |4] we show the radial distribution of the thermal 
Bremsstrahlung X-ray flux for a SNR expanding in a medium 
with density n = 0.1 cm^'^. The SNR age is 1000 yr and the ex- 
plosion energy is Esn = 10^^ erg. We neglect the effect of CR 
acceleration at the SNR shock (i.e. w = 0). The solid, dashed 
and dot-dashed lines refers to the emission profiles at 0.5, 1 and 
2 keV, respectively. The emission is roughly constant in the cen- 
tral region and exhibit a pronounced peak close to the position 
of the shock, indicating the presence of a shell-type morphology 
in thermal X-rays. 



3. Results 

Here we solve numerically Eq. |4]and|5]to obtain the electron 
and proton temperature distributions inside a SNR. The depen- 
dence of the projected temperatures Te.p{x — 0.7) (hereafter 
referred to as simply Tg ,,) over parameters such as the upstream 
gas density n and the supernova explosion energy Esn can be 
described as power laws for a broad range of parameter values 
in = Om..lcm-^,t^ 125.. 2000 yi, Esn = lO'^.-S-lO^i erg, 
w — 0..0.8). Parameterizations for the electron temperature are 
given in tables IaTI and IaTSI (see Appendix [All. Here we briefly 
summarize the main characteristics of the solution. 

Power law approximations provide a good fit to the exact 
solution except for the cases in which T^ ^ Tp. Fig. |5] illus- 
trates this situation: the electron and proton temperatures are 
shown as a solid red and dashed green curve respectively, for 
the two cases Tg ^ Tp (left panel) and Tg ^ Tp (right panel). 
The left panel shows the case of SNR of the age t = 2000 yr 
with significant CR pressure (w = 0.8). The explosion en- 
ergy is Esn — lO^^erg. Under these assumptions for densi- 
ties n ;» 0.1 cm^"^, we get almost equal electron and proton 
temperatures. The dependence of these temperatures from the 
density deviates from a power law scaling. The right panel of 
Fig. prefers to a SNR with age t — 500 yr, total explosion en- 



ergy E = lO^^erg and negligible acceleration of particles at the 
shock {w = 0). In this case T^ ^ Tp and as it is evident from 
the figure temperatures scale as a power laws of the density. 
In general, results can be summarized as follows: 

- from t = 125 yr to i = 2000 yr and for w < 0.5 the 
dependence of T^ on gas density n can be approximated 
with a single power law Tg oc n^ with index in the range 
b = 0.15 - 0.23 for n=0.01..1cm-3. 

For t > 1000 yr and w > 0.5, the electron temperature Tg, 
in general, can not be approximated with a power law of the 
gas density (see AppendixlAlfor details). 

- The dependence of T^ on the total explosion energy Esn 
can be approximated as a power law Tg oc £'5^ with index 
in the range b = 0.16 — 0.31. This is valid for the whole 
considered range of parameters (see Appendix lAli. 

It has been argued that some additional mechanism of elec- 
tron heating based on the generation of plasma waves must 
operate in SNR s (Caraifl & PacadoDoulos 1988; Laming 200ll 
iGhavamian eTall KOO?: Rakowski et al. 2008). 

Apriori, two possibilities for this mechanisms could be con- 
sidered. They might operate only in the vicinity of the shock 
(■Cai-gill & PaDadoi30ulos"l988'; 'Laming"200lt IGhavamian et al.l 
2007; Rakowski et al. 2008), either, similar to Coulomb heating 
mechanism, operate all way downstream. Some physical justifi- 
cation f or the second type mecha nis ms could be found in|Lamina 
(I2004t) ; iLaming & Lepril (l2007h and lSharma et al.1 (l2007l) . were 
analogous processes for the solar wind and hot accretion flows 
were considered. 

Formally speaking, mechanisms operating very close to the 
shock, implies changing of boundary conditions for our equa- 
tions. We will consider them in details for SN 1006 in the sec- 
tion |4] Results for mechanism that operate all way downstream 
are briefly summarized below. 

If we assume for simplicity that this mechanism has the same 
temperature dependence as the Coulomb one, but is k times 
faster, then the heating of electrons can be described by an equa- 
tion similar to Eq. |4l 



dt 



k^ 



n 



7 



eq 



n dt 



By solving this equation we found that the electron temperature 
scales with k as T^^k) ~ fc0-8-0-9. 

A comparison between the predictions of our model, which 
considers only Coulomb heating, and the observed SNR electron 
temperatures can serve as a tool to estimate whether such addi- 
tional heating mechanisms are needed at all, and, if so, at which 
level. 



4. Application of the model to the supernova 
remnant SN 1006 

In order to compare the predictions of our model with obser- 
vations we consider here the well known SNR SN 1006. The 
choice of this object is justified by the fact that it is a very well 
studied SNR, for which a relatively accurate determination of 
parameters such as distance, shock velocity and age, is available. 
SN 1006 is an historical SNR with age w 1000 yr and be- 
lieved to be in the Sedov phase of its evolutio n. After reviewing 
the preexisting literature. lDubner et al.l (1200 2') concluded that the 
distance to SN 1006 lies in the range 1.4-2 kpc. More rec ently, a 
slightl y larger distance of 2.18 kpc was inferred by Winkl er et al.l 
(l2003h . who compared the measured optical proper motion with 
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Fig. 5. The left panel shows the electron and proton temperatures for a SNR of age t = 2000 yr, CR acceleration efficiency 
w = 0.8 and explosion energy Esn = lO^^erg. The right panel shows analogous dependencies for a SNR for t ~ 500yrs, w — 0, 
EsN = lO^^erg. In both cases projected temperatures are plotted, for a distance from the SNR centre equal to r/r shock =0.7. 



the expansion velocity of 2890 ± 100 km/s determined from Ha 
lines observations. This velocity is i n a good agreeme nt with 
analogous velocity measurements bv lSmith et al.l (11991') which 
give V = 220 - 3500 km/s. At th is distance the SNR appar- 
ent size 15' (jRothenflug et al.ll2004l) translates into a shock ra- 
dius equal to rghock = 9.6 pc. A value for the ambient gas 
density equal to n ~ 0.15 — 0.25 cm~^ ( i ri the n orth-western 
region) has been obtained by 'Acero et aP (l2007h from XMM 
X-ray observations. A slightly higher values for the density 
n — 0.25 — 0.4 cm~'^ was obtai ned from Ha lines observations 
with HST (iR avmond et ani2007h . 

By assuming that the evolution of SN 1006 can be de- 
scribed by the Chevalier self-similar model described in Sec. |2] 
it is possible to derive a supernova explosion energy equal to 
EsN~2-K 10^1 erg. 

Results from our model for SN 1006 are shown in Fig. |6] 
where the electron and proton temperatures behind the SNR 
shock wave are plotted as a function of the remnant age. For 
the actual age of the remnant, t k, 1000 yr, ambient gas density 
n = 0.2 cm~'^, explosion energy E = 2 ■ 10^^ erg and w — 0.1, 
our model gives a value of the proton temperature equal to Tp ^ 
17 keV, and an electron temperature equal to T^ = . 6 keV . 
The observational values are adoptecQ from IVink et all (l2003h 
for the proton and electron temperatures are also indicated in 
Fig. |6] together with the error bars (3ct for the proton temper- 
ature and 90% confidence level for the electron temperature). 
A bit lower prot on temperature Tp = 1.8 • lO^K « IGkeV 
was obtained by iKorreck et al.l (|2004|) . It is clear from Fig. |6] 
that our model predicts a proton temperature which is in agree- 
ment with observations, while the prediction for the electron 
temperature falls a factor of 2-3 below the observational value 
of Te « 1.5 ± 0.2 keV (IVink et all l200l . (XMM). A simi- 
lar value for the electron temperature of Te = 1 . 5 — 1 .7 keV 
has been independently measured by lAcero et aP (l2007l) with 
XMM. Nevertheless there is some doubt in the electron temper- 
at ure measurements. Lower (0.6-0.7keV) values were obtained 
bv lLongetal] (l2003h with CHANDRA for NW-1 region, while 
other NW regions appear to be ejecta dominated. Since XMM 



' Note, that in IVink et all J2003h not the proton, but oxygen tem- 
perature is measured directly. We adopt here Tp — mp/moTo ~ 
l/WTo «33keV 



has worse spatial resolution than CHANDRA, ejecta dominated 
regions can contribute significantly to its spectra. Moreover, in 
I Acero et aP (l2007h it is shown that including synchrotron emis- 
sion component can reduce observed electron temperature. 

If density and velocity vary in the range n = 0.15 — 
0.5 cm^'^, V = 2400 — 3500 km/s, respectively, we estimated the 
corresponding range of electron temperatures predictions from 
the model to be Tg = 0.5 — 0.9 keV. This value s are in the agree- 
ment with ones from the model in IVink et al.l (2003), see their 
Fig. 2 . The shaded regions in Fig. |6] represent the predictions 
from our model by taking into account the above mentioned un- 
certainties on SNR parameters. 

Thus, the predicted ele ctron t e mpera ture agrees well with 
fits to N W- 1 region in Long et al. (l2003l). but smaller, than th e 
temperatures obtained bv lVink et al.l(l2003h : lAcero et alj(l2007h . 
Part of the discrepancy between the measured and expected elec- 
tron temperature might be explained by the uncertainties in the 
determination of SNR parameters such as the external gas den- 
sity or the shock velocity. 

Dependence of the electron and proton temperature on accel- 
eration efficiency for SN 1006 is shown in Fig.|7] Both tempera- 
tures quickly decreasing with increasing w. The result shown in 
Fig.|6]corresponds to w = 0.05. 

Recent detection of S N 1006 NE and SW re gions at TeV 
gamma rays by H.E.S.S. dAharonian et al.l l2010h implies that 
CRs protons and/or electrons are accelerated at the SNR shock. 
The density rio « 0.085 cm"^ ( Katsuda et al. 2009 ) in this re- 
gions is lower than in NW and thus they can be in early Sedov 
or late free-expansion phase. The measurements of expansion 
index by Katsuda et al. (2009) give m — 0.54 ± 0.05, i.e. con- 
siderably less than if it were in a free expansion (m — 1), but 
higher than for an adiabatically expanding SNR {rn — 0.4). For 
an estimation of the electron temperature we consider this SNR 
to be in Sedov phase. Under the assumption that the observed 
TeV emission has an hadronic origin it is possible to estimate 
the CR acceleration effi ciency to be roughly equal to w w 0.2 
dAharonian et al.ll2010l) . For efficiency range w — 0.1. .0.2 our 
model gives electron temperature Tg « 0.47 — 0.52 keV. 

The expected electron temperature can be further increased 
if additional heating mechanisms for electrons operate at SNR 
shocks. It has been suggested that these mechanisms might be re- 
lated to the presence of plasma waves upstream of the shock that 
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Fig. 6. Solid and dashed lines show the electron and proton tem- 
perature for SN 1006 as a function of the SNR age. Shaded re- 
gions represent the model uncertainty (see text for details). Data 
points refer to the observed va l ues of Tg and Tp (dots with error 
bars) adopted from lVink et al.l (l2003h . see text for the details. 



can effectively heat electrons d Cargill & Papador)oulos 1988^, 
iLamind 120011: iGhavamian etaLl ,2007.; .Rakowski et al., ,2008,) . 
The fact that such heating is expected to happen upstream of 
the shock is suggested by optical observations of Ha lines from 
a number of SNR shocks (see Ghavamian et al. 2007, and refer- 
ences therein). These observations show that the electron tem- 
perature might be akeady significantly heated above the ex- 
pected value of Te — {me/mp) Tp immediately after the shock, 
since the width of Ha filaments is of the order of « 10^'' cm. It is 
important to stress that, on the other hand, the temperatures mea- 
sured from X-ray observations of SNRs refers to a quite broad 
region downstream of the shock (e.g. ^10%rs/i for SN 1006). 

In order to estimate the effect of additional heating hap- 
pening upstream of the shock, we repeated our calculations for 
SN 1006 by assuming that the ratio of the proton to the electron 
temperature immediately after the shock is equal to T^/Tp = 
1/15, instead of the canonical value T^/Tp — me/mp w 
1/1800 assumed in Fig. |6] Such value is consistent with the 
range of values derived from H a measurements T^jTp < 0.07 
(IGhavamian et al.l l2002l 12007). Thus, this is equivalent to as- 
sume that the additional heating increases the electron tempera- 
ture at the shock by a factor of w 100. For such a choice of the 
boundary condition our model predicts a value for the effective 
electron temperature (the one that is measured from X-ray ob- 
servations) equal to 1 .4 keV. This value is higher than the one 
obtained by the minimal heat ing model shown in Fig. |6] and in 
agreement with observations (IVinket al.ll2003l) . It is important 
to stress that an increase of about two orders of magnitude of the 
heating of electrons at the shock (which leads to T^/Tp = 1/15 
immediately downstream of the shock), is reflected by a moder- 
ate increase (a factor of w 2 only) of the predicted effective tem- 
perature which can be measured from X-ray observations. This 
is due to the fact that the electron temperature well downstream 
of the shock is mainly determined by an equilibrium between 
Coulomb heating and adiabatic losses and depends only mildly 
on the boundary condition for T^/Tp at the shock. 



Electron temperature 
Proton temperature 



0.4 0.5 0.6 0.7 

Acceleration efficiency 



Fig. 7. Solid and dashed lines show the electron and proton tem- 
perature for SN 1006 as a function of acceleration efficiency w. 



5. Discussion and conclusions 

In this paper, we studied the evolution in time of the proton 
and electron temperatures behind the shock waves of SNRs in 
the Sedov phase of their evolution. To descri be the evolution o f 
the SNR we adopted the model described by Chev ahen (Il983h . 
which generalizes the Sedov solution by allowing the accelera- 
tion of CRs at the SNR shock. Immediately after being shocked, 
the proton gas is heated to the temperature given by the mod- 
ified (by the inclusion of CR acceleration) Rankine-Hugoniot 
relations, while the electron temperature was assumed to be a 
factor of mp/m,e lower. Downstream of the shock electrons can 
be further heated due to Coulomb collisions with the hot protons. 

Since Coulomb heating is unavoidable and is the only heat- 
ing process included in our calculations, our results has to be 
intended as estimates of the minimal electron temperatures ex- 
pected in SNRs. 

We applied the model to the case of SN 1006. The model pre- 
dicts a value of the electron temperature equal to T^ — 0.6 keV, 
roughly a factor of 2 smaller than the value measured from X-ray 
data, which is Tg — 1.5 keV. The discrepancy may be reduced 
by taking into account the uncertainties in the determination of 
the SNR parameters such as the shock speed or the external gas 
density. 

In fact, other heating process based for example on the gen- 
eration o f plasma waves upstrea m of the shock may play a role 
(see e.g. iRakowski et al.l l2008l) . Such additional collisionless 
processes have been sometimes invoked to explain the observed 
electron temperatures, since Coulomb hea ting alone seemed to 
be too slow in certain cases (lLaming| l2000). In this respect, our 
model can serve as a useful tool to estimate the order of mag- 
nitude of such extra heating, when that is required. For the con- 
sidered case of SN 1006, we showed that the discrepancy of a 
factor of « 2 between the expected electron temperature and the 
one measured from X-ray observations can be indeed compen- 
sated if an additional heating mechanism operates upstream of 
the shock. This mechanism has to be very effective and able 
to increase the electron temperature immediately downstream 
of the shock (e.g. before Coulomb heating starts to operate) of 
about two orders of magnitudes above the canonical value of 
Te = (me/mp) Tp. 
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Appendix A: Power law fits for tlie eiectron 
temperature 

In this Appendix we present the coefficients for the power law 
fits of the effective electron temperature T^ as defined in Sec. 12. 31 
For fixed values of the SNR age t and CR acceleration efficiency 
w, the electron temperature T^ is fitted as a power law function 
of the ambient gas density n and explosion energy Esn- Fits 
parameters are given for values of w in the range 0..0.8 and for 
vlaues of the SNR age in the range 125.. 2000 yrs. 

Table lA.ll gives the values of the parameters for the power 
law fits: 

Te(n)/lkeV = a ■ (n/lcm-^)''. 

All fits have an accuracy better than «10% for values of the am- 
bient density in the range 0.01..1 cm^'^. For some sets of param- 
eters a power law fit could not give a satisfactory approximation 
(better than 10%) for Te- For these cases the electron tempera- 
ture Te{n) is plotted in Fig. A.l and Fig. A. 2. Each curve in the 
figures represents Te{n) for different values of explosion energy, 
see figures captions for the details. 

Table IA.2I gives the values of the parameters for the power 
law fits: 

Te(£;)/lkeV = a ■ (£:/10^^erg)^ 

The accuracy of fits is better than « 10 % when the value of the 
ambient density is in the range n — 0.01..1 cm"-^ and the CR 
acceleration efficiency is int he range w = 0..0.8. 



D. Malyshev et al.: A simple model for electron plasma heating in supernova remnants. (RN) 



E^ 


^QbU 


3- 


^QbO 


5- 


10"" 


7- 


^QbU 


10^^ 


3- 


10'^^ 1 


w=0 1 


4 


a 


b 


a 


b 


a 


b 


a 


b 


a 


b 


a 


b 


125 


0.75 


0.217 


0.89 


0.198 


0.96 


0.194 


1.01 


0.187 


1.06 


0.178 


1.25 


0.145 


250 


0.70 


0.230 


0.84 


0.228 


0.92 


0.227 


0.97 


0.225 


1.03 


0.224 


1.23 


0.216 


500 


0.65 


0.231 


0.78 


0.234 


0.86 


0.230 


0.91 


0.232 


0.96 


0.230 


1.16 


0.231 


1000 


0.59 


0.215 


0.72 


0.221 


0.79 


0.226 


0.83 


0.224 


0.89 


0.225 


1.07 


0.227 


2000 


0.50 


0.187 


0.63 


0.199 


0.69 


0.204 


0.74 


0.207 


0.79 


0.209 


0.97 


0.216 


w=0.3 1 


125 


0.57 


0.229 


0.68 


0.225 


0.74 


0.222 


0.79 


0.220 


0.83 


0.217 


0.99 


0.203 


250 


0.44 


0.226 


0.53 


0.228 


0.57 


0.229 


0.61 


0.230 


0.64 


0.230 


0.94 


0.228 


500 


0.49 


0.220 


0.59 


0.224 


0.65 


0.226 


0.68 


0.227 


0.73 


0.227 


0.88 


0.229 


1000 


0.42 


0.198 


0.52 


0.208 


0.58 


0.211 


0.62 


0.213 


0.66 


0.215 


0.80 


0.221 


2000 


0.33 


0.148 


0.43 


0.170 


0.48 


0.178 


0.52 


0.183 


0.56 


0.187 


0.70 


0.199 


w=0.5 1 


125 


0.46 


0.231 


0.56 


0.230 


0.61 


0.228 


0.65 


0.228 


0.69 


0.226 


0.82 


0.221 


250 


0.44 


0.226 


0.53 


0.228 


0.57 


0.229 


0.61 


0.230 


0.64 


0.230 


0.77 


0.231 


500 


0.39 


0.212 


0.48 


0.218 


0.52 


0.220 


0.56 


0.222 


0.59 


0.223 


0.72 


0.226 


1000 


0.32 


0.180 


0.41 


0.193 


0.45 


0.199 


0.49 


0.202 


0.52 


0.205 


0.64 


0.213 


2000 


- 


- 


- 


~ 


- 


~ 


- 


~ 


- 


- 


0.55 


0.178 


w=0.8 1 


125 


0.30 


0.226 


0.36 


0.228 


0.39 


Q.ll^ 


0.41 


0.230 


0.44 


0.230 


0.53 


0.231 


250 


0.27 


0.213 


0.33 


0.219 


0.36 


0.221 


0.38 


0.222 


0.41 


0.223 


0.49 


0.227 


500 


0.23 


0.183 


0.28 


0.196 


0.31 


0.201 


0.34 


0.204 


0.36 


0.207 


0.44 


0.214 


1000 


- 


- 


- 


- 


- 


- 


- 


- 


- 


- 


0.38 


0181 


2000 


- 


- 


- 


- 


- 


- 


- 


- 


- 


- 


- 


- 



Table A.l. Fit for the electron temperature Te for SNR inner regions (i.e. at projected distance r jr shock ^ 0.8) as function of outer 
density Te(n)/ Ike V ~ a(n/lcm^'')'', for different explosion energies, if applicable. SeeFig.A.l,Fig.A.2 otherwise. Uncertainties 
of all fits are less than «10% 
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Fig. A.l. Electron temperature Tg for SNR inner regions (i.e. at projected distance r/r shock ^ 0.8) for parameters for which 
powerlaw fit is not applicable (see table lA. il l. Left panel: TQ{n) dependence for t=2000, w=0.5. Each curve represent different 
explosion energy. Right panel: Te(n) dependence for t=1000, w=0.8. Each curve represent different explosion energy. 
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Table A.2. Fit for the electron temperature Te for SNR inner regions (i.e. at projected distance r jr shock ^ 0.8) as function of 
explosion energy T^ (i?)/lkeV = a ■ (i?/10'^^erg)^ for different outer gas densities. Uncertainties of all fits are less than 10% 
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Fig.A.2. Same as Fig.A.l fort=2000, w=0.8 
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